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Интерпретация динамики взаимодействия 
судна с экстремальной волной в бортовой 
интеллектуальной системе 


Обсуждается проблема контроля критической ситуации, связанной с воздействием на судно экстремальной 
волны. Особое внимание уделяется анализу особенностей взаимодействия в условиях неопределенности и 
неполноты исходной информации. 


Введение 

Создание интеллектуальных систем (ИС) новых поколений в значительной 
степени определяется уровнем интеллектуальности, достижение которого открывает 
возможности проектирования систем, функционирующих в условиях неопределен- 
ности [1-27]. Особенно важно это в нештатных ситуациях, которые возникают при 
воздействии экстремальных волн огромной разрушительной силы. При взаимодей- 
ствии с внешней средой судно, как сложный динамический объект (ДО), подвергается 
возмущениям, которые заранее трудно предусмотреть при создании системы. Это 
накладывает дополнительные требования к обработке измерительной информации. 
Дефицит времени на принятие решения по обеспечению безопасности ДО приводит 
к необходимости выделения только той информации, которая требуется для выпол- 
нения основной задачи управления и принятия решений. Возникает проблема пред- 
варительного отбора и анализа информации, необходимой для реализации меха- 
низма логического вывода и выработки практических рекомендаций. Решение этой 
проблемы требует создания эффективных алгоритмов и процедур обработки инфор- 
мации, поступающей с датчиков измерительной системы. Рассматриваемая ИС обеспе- 
чивает оперативный контроль текущей ситуации. На основе данных динамических из- 
мерений необходимо построить эффективный прогноз и реализовать управление. Вся 
другая информация имеет смысл только на этапах предварительного анализа при по- 
строении ИС [3], [5]. 

Ниже рассмотрена одна из наиболее сложных и трудноформализуемых задач, 
решаемых при интеллектуальной поддержке оператора ИС - оценка поведения судна 
при воздействии экстремальной волны в сложной динамической среде, характеризу- 
ющейся полной неопределенностью как объекта моделирования, так и динамики 
внешней среды. Крупные волны в Мировом океане встречаются достаточно часто и 
представляют грозную опасность для судоходства. Вместе с тем среди экстремаль- 
ных волн можно выделить гигантские волны огромной разрушительной силы, пара- 
метры которых не соответствуют общепринятым представлениям о ветровых волнах. Это 
так называемые «необычные» волны или «волны-убийцы». В англоязычных публика- 
циях такие волны получили название «НеаК \’ауез», «гогие \ауез» или даже «КШег- 
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У’ауез», т.е. «волны-убийцы». Под «волной-убийцей» понимают внезапно возникаю- 
щую интенсивную волну, которая значительно превышает (в 2 и более раз) значи- 
тельную высоту фонового волнения (рис. 1). 

Внезапность и большая энергия таких волн неоднократно приводили к гибели 
судов и морских сооружений. Долгое время считавшиеся фантазией мореплавателей, 
«волны-убийы» стали объектом серьезного изучения после участившихся случаев 
регистрации подобных волн на обширных акваториях океанов и морей. Наиболее 
известным примером «волны-убийцы» стала так называемая «новогодняя волна», 
зарегистрированная на нефтяной платформе Экофиск (56,5 с.ш., 30,2 в.д.) в Север- 
ном море 1 января 1995 года. Эта необычная волна имела высоту 25,6 м, гребень 18,5 м, 
ложбину 7,1 м. 


высокопроизводительные вычисления и моделирование в сложных динамических средах 
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Рисунок 1 — Параметризация необычной волны по данным измерений. 
(Черное море, 16 декабря 2000 г.) (а); 
динамическая картина взаимодействия в экстремальной ситуации (6) 


В течение длительного времени наиболее распространенной точкой зрения яв- 
лялась некоторая исключительность акваторий, где возможны необычные волны, 
образующиеся, главным образом из-за наличия встречного течения для волноопас- 
ных направлений, резкого свала глубин и т.п. Однако по мере накопления данных о 
встречах с «волнами-убийцами» географический фактор теряет свою уникальность и 
уступает место принципу возможности встречи с этой волной на всем пространстве 
Мирового океана. 

Исследование воздействия «волны-убийцы» на суда и плавучие технические 
средства освоения океана требует кардинальной модификации математической мо- 
дели за счет включения в нее механизмов описания дополнительных эффектов, не 
учитываемых обычной гидродинамической моделью взаимодействия в условиях 
существенной неопределенности задачи. Физические закономерности и необычная 
структура «волны-убийцы» порождают ряд проблем, решение которых возможно 
только на базе новых концепций и упрощающих предположений. 

Формулировка теоретической базы исследования динамики взаимодействия с 
«волной-убийцей» ведется с учетом физических картин взаимодействия судна с внеш- 
ней средой. Это позволяет прояснить особенности поведения судна и избежать мно- 
гих трудностей при решении проблемы методологического барьера — нелинейности, 
многомерности, многосвязности [3], [9], [11], [13], [17]. На основе такой формализа- 
ции может быть построено математическое описание поведения судна как нелиней- 
ной динамической системы, определяющей процесс восприятия судном возникаю- 
щих динамических нагрузок. 
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Нечеткая формальная система 


Формализация понятий и знаний о динамике взаимодействия судна с внешней 
средой при реализации системы информационной поддержки оператора бортовой 
ИС достигается на основе теории нечетких множеств [1]. Нечеткие отношения 

К=А-> В, (1) 
заданные между входной областью Х и выходной областью У в виде нечеткого под- 
множества прямого произведения ХхУ, определяются по формуле 


в=>.>.{ (ль), (2) 


1 = 
где Х = {х.....х,} — область посылок (антецедент импликации); У = {у1,...,Уш — 
область заключений (консеквент импликации); ив(хьу;) - функция принадлежности 
(ФП) (хьу} нечеткому отношению К; 4, (х,, у,) [0,1]; знак У означает совокупность 
(объединение) множеств. 


Для продукционных правил, использующих нечеткие множества А(АСХ) и 
В(ВсУ), один из распространенных способов построения нечеткого соответствия 


имеет вид: 
в=Ахв=У У (инк) ль, ,)}} 
у (3) 
ИЛИ 
мабсу) = наС®) ^ ив(у) = шина ^ нв(У), (4) 


где ид(х) и ив(у) — ФП элементов х,у множествам А и В. 

Алгоритмы нечеткого вывода в основном различаются видом используемой 
операции нечеткой импликации. В бортовом вычислительном комплексе находят при- 
менение различные модели нечеткого вывода, основанные на расширении принципа 
композиции Л. Заде [27]. В основном это модели Мамдани и Сугено [3]. 

Рациональный учет особенностей поведения судна под воздействием внешних 
возмущений требует разработки системы критериальных оценок, обеспечивающих 
проверку безопасности судна в рассматриваемой экстремальной ситуации в виде 
глобальной и локальной систем нормирования [5]. Глобальная система критериев 
включает национальные и международные требования, которые обеспечиваются не- 
зависимо от типа судна. Локальная система разрабатывается в процессе создания си- 
стемы знаний и учитывает характерные особенности исследуемого судна, основан- 
ные на более совершенных методах формализации информации с учетом ее неполно- 
ты и неопределенности в рамках концепции «мягких вычислений» [27]. Однако 
контроль динамики судна при воздействии волны-убийцы имеет свои особенности, а 
обычные критерии остойчивости не могут гарантировать безопасность судна в этой 
экстремальной ситуации. 

Следует отметить, что ситуация, связанная с ударом волны-убийцы, имеет неко- 
торое сходство с ситуацией возникновения «брочинга» [3], [12], [21] в том смысле, 
что никакие действия судоводителя в момент возникновения этой ситуации не могут 
привести к обеспечению безопасности. Поэтому при формализации знаний в борто- 
вых ИС целесообразно предусмотреть логическую систему, обеспечивающую прог- 
ноз возникновения волны-убийцы в заданных внешних условиях, а также сценарии 
поведения судна под ее воздействием. 

Рассмотрим некоторые общие соображения об особенностях преобразования 
информации при построении нечеткой формальной системы, обеспечивающей при- 
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нятие решений на основе анализа альтернатив. При построении формальной системы 
на основе имеющейся информации введем в рассмотрение множество векторов, 
характеризующих конкретную альтернативу (ситуацию): 
В": К = (жь,..., хи) > В", (5) 

гдехи,..., хи — признаки, выделенные в процессе анализа измерительной информации. 

Таким образом исходная информация о динамике ветроволновых полей может 
быть представлена структурой, определяемой пространство принятия решений. Число 
возможных значений лингвистических переменных (мощность терм-множеств) опре- 
деляется в процессе предварительного анализа при построении матрицы стратеги- 
ческих решений [3], [23]. Формализованное описание лингвистической переменной 
(показателя состояния динамической системы) представляет собой упорядоченную 
совокупность: 

Вы = (т, у, ших}, Х >, (6) 
где Аз — наименование состояния; Т — терм-множество (множество состояний) лингви- 
стической переменной: \ — имя нечеткой переменной; у и ш(х)/х — ее базовая шкала и 
функция принадлежности; Х — базовая шкала лингвистической переменной. 

Формализация пространства принятия решений как нечеткой системы осущест- 
вляется на базе представления знаний об исследуемой предметной области в виде 
совокупности нечетких формальных систем и переходов между ними [23]. Нечеткая 
формальная система учитывает динамику изменения предметной области и пред- 
ставляется в виде матрицы переходов, строками которой являются пары 


(ИЗь <Ть $25 1=1,.. 3] = 1..5), (7) 


где и91 — расширенная нечеткая формальная система; <Т:,3}> — пара, задающая мак- 
симинную операцию, используемую в основе композиционного правила вывода 
(обобщенный тойих ропеп5); Т,— оператор минимума (Т-норма); 5] — оператор мак- 
симума (Т-конорма). 

Условия переходов определяются при разработке формальной системы. При 
этом возможны переходы от одной расширенной нечеткой формальной системы к 
нескольким подобным системам, которые могут создаваться и удаляться в процессе 
моделирования динамики системы, определяющей конкретный ветроволновой режим. 


Интерпретация поведения судна 
в рассматриваемой экстремальной ситуации 


Предлагаемый ниже подход к формализации задачи основан на использовании 
синергетической парадигмы [6], [8], [10], [19], [20], [24]. Интерпретация поведения 
судна, как сложной динамической системы (рис. 2), проведена в терминах синер- 
гетической теории управления [8]: 

1. Предполагается, что судно, как динамическая система, под воздействием 
внешнего возмущения «стартует» из некоторого исходного состояния, характеризуе- 
мого множеством начальных условий фиксированной размерности и через опре- 
деленное время неизбежно попадает на формирующийся аттрактор — притягивающее 
множество, которое имеет значительно меньшую размерность фазового пространства. 

2. Эволюция сложной динамики — процесс сжатия фазового пространства, отра- 
жает особенности самоорганизации системы, развитие которой происходит в течение 
чрезвычайно ограниченного периода времени и может быть описано модифициро- 
ванными математическими моделями поведения судна в различных фазах взаимо- 
действия. 
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3. Расчетная схема использования синергетической парадигмы при оценке 
поведения судна при воздействии экстремальных волн может быть построена на базе 
стандартных соотношений, используемых в синергетической теории управления, и 
математических моделей, описывающих различные фазы эволюции системы в 
процессе движения к целевому аттрактору. 


Результаты вычислений 


Процедура. основанная на подходе Процедура, основанная на подходе 
С.Н.Благовещенского и А.Н.Холодилина Н.Б.Севастьянова 


Качка судна Компоненты математической модели 
На больших глубинах: 
О [пей = 0,0328 рад^2 О _\И нев = 0,0874 рад^2/с^2 
Треа_0.1% = 38,4 град УТреа_0.1% = 1,094 рад 
На мелководье 
0_фаа_пу = 0,0158 рад^2 О_\ ева пм = 0,0442 рад^2/с^2 
Треа_0.1% = 26,6 град УТневз_0.1% = 0,778 рад 


Момент инерции судна [1х+144) = 14,4 тс^2/м 

Козффициент демпфирования = 016 1/с 

Максимальное значение кренящего момента = 38,66 тм 

Кренящий момент начинает действовать в момент ( = 3,8 с, 

когда счано проходит прямое положение при двихж в напр бега волн 


Силовое воздействие разрушающейся волны Изменение угла крена 


Сила удара разруш волны = 25,47 т 

Время действия цдара = 1,24 с 

Импульс удара = 15,74 тс 

Аппликата центра цдара волны = 2,3 м 

Угл скорость качки после удара = 1/11 с^[1] 

00,2 0,54 0,94 1,33 1,73 212 2,51 2,91 3,3 38494 4,34 473 


Кинетическая энергия 
Зер, с 


Кинетич энергия судна после удара = 17.52 тм —— Наклонение судна в условиях разруш волнения 
Потенц энергия судна = 1919 тм —— Наклонение судна без воздействия разруш волн 


Судно выдерживает удар разрушающейся волны. Внимание! Судно НЕ выдерживает чдара разруш волны! 


Вспомогательные данные | : Закрыть 


Рисунок 2 — Результаты моделирования при времени действия удара 
разрушающейся волны (время действия удара 1,24 с) [3] 


Особенность использования синергетического подхода при анализе динамики 
судна в критической ситуации, связанной с воздействием экстремальных волн, со- 
стоит в универсальности методов исследования. Благодаря бифуркации нелинейность 
процессов самоорганизации определяет многовариантность путей развития системы. 
При этом принцип подчинения сужает диапазон величин, изменение которых с тече- 
нием времени определяет процесс самоорганизации. Несмотря на то, что исследуемая 
динамическая система имеет сложный характер и требует учета многих переменных, 
применение принципа самоорганизации позволяет ограничиться лишь небольшим 
числом самых важных параметров — параметров порядка, чтобы эффективно наблю- 
дать за развитием системы. 


Геометрическая и аналитическая интерпретация 
ситуации с позиций теории катастроф 


Общая расчетная схема оценки динамики судна в рассматриваемой задаче на 
основании сформулированной гипотезы представляет собой достаточно громоздкую 
расчетную схему, реализуемую в мультипроцессорной вычислительной среде и пре- 
дусматривающую последовательную проверку безопасности судна в трех последова- 
тельных стадиях развития экстремальной ситуации. Эволюция системы в этих ста- 
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диях характеризуется постепенным уменьшением числа степеней свободы и форми- 
рованием аттракторных множеств. В завершающей стадии развития эволюции при 
реализации эффективного управления может быть достигнуто возникновение целе- 
вого аттрактора. Если динамические нагрузки от удара экстремальной волны окажутся 
катастрофическими, то в процессе эволюции будет реализована критическая ситуация, 
приводящая к опрокидыванию судна. 

В качестве иллюстрации динамики судна, воспринимающего удар экстремальной 
волны, рассмотрим моделирование этой ситуации на основе модифицированной мо- 
дели, разработанной на основе методов теории катастроф [2], [4], [7], [18], [22]. Дина- 
мическая модель катастрофы, отображающая поведение сложного объекта, может быть 
представлена в виде кортежа: 


ОМ(Са® = < Н(Са®, А(Сай >, (8) 


где Н(Саб) — компонента, интерпретирующая динамическую среду на основе геоме- 
трической модели катастроф; А(Саб) — компонента, интерпретирующая динамиче- 
скую среду на основе аналитической модели катастроф. 

Геометрическая интерпретация физической картины взаимодействия на основе 
компоненты Н(Сай) имеет вид: 


Н(Са® = < В(6, ®), 97((6, ®, С(6, ®) >, (9) 


где В(0,1) — бифуркационное множество, представляющее собой образ особого мно- 
жества в пространстве управляющих параметров; ©7(0,1) — множество, отображающее 
динамическую среду, интегрирующую влияние внешних возмущений и особен- 
ностей динамики объекта; С(9,1) — множество, отображающее структурные измене- 
ния в форме подводной части корпуса, вызванные непрерывным изменением конфи- 
гурации действующей ватерлинии. 

Аналитическая интерпретация физической картины взаимодействия на основе 
компоненты А(Сай) имеет вид: 


А(Сай) = < А(СВ), А((РН), АС\Е) >, (10) 


где А(СК) — математическая модель, отображающая оценку динамической среды на 
основе критериальных соотношений, разработанных в рамках стандартных алгоритмов 
исследуемой проблемной области; А((РН) — математическая модель, отображающая 
оценку динамической среды на фазовой плоскости с помощью теории стохастических 
систем, теории детерминированного хаоса и синергетической парадигмы; А(МЕ) — 
математическая модель, отображающая оценку динамической среды на основе нейро- 
нечеткой системы. 

Геометрическая интерпретация сложной динамической системы для модифи- 
цированной модели сборки представлена на рис. 3. 

Применим типичную для рассматриваемого случая взаимодействия катастрофу 
сборки и рассмотрим возможные сценарии развития ситуации в рамках синергети- 
ческой парадигмы: 

— удар экстремальной волны при полном погружении судна, которое можно 
рассматривать как некоторое подводное судно с необычной формой корпуса и исполь- 
зовать модель взаимодействия, изображенную на рис. 4 а; 

— развитие стремительного дрейфа вследствие удара экстремальной волны (рис. 4 6); 
результирующая нагрузка от гидродинамических сил в этой ситуации образует кре- 
нящий момент относительно ЦТ системы; величина этого момента непрерывно уве- 
личивается в процессе развития дрейфа, что приводит к возникновению катастрофи- 
ческого крена или опрокидыванию судна. 
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В(6,‚1) - бифуркационное множество, 
представляющее собой образ особого 
множества в пространстве управляющих 
переменных 


С2(60,;1) - многообразие, определяющее 
динамическую среду, интегрирующую 
влияние внешних возмущений и 
особенностей динамики объекта 


С(0,4) - множество, отображающее 
структурные изменения в форме 
подводной части корпуса 


Рисунок 4 — Топологические картины взаимодействия судна с экстремальной 
волной: а) — в момент удара; 6) — после удара 


Как видно рис. 4 а), при таком сценарии происходит полная перестройка топо- 
логической картины катастрофы сборки. Характерные точки С, Си, по, Р и область 
бифуркационного множества В изменяют свое положение. При этом восстанавлива- 
ющий момент создается за счет положительной метацентрической высоты М(0) = 
= Ороз110. В идеальном случае (судно с круговыми обводами) метацентр то и центр 
величины (ЦВ) Со совпадают и находятся в центре круга. Если судно выдерживает 
удар, то возникает ситуация, описанная в следующем сценарии. 

После восприятия удара снова происходит перестройка топологической картины 
катастрофы сборки (рис. 4 6)). ЦТ системы располагается в точке Со*, и множество, 
отображающее динамическую среду, приобретает вид С(70(0)* (светлая область). 
Диаграмма остойчивости М(0) в этих условиях формируется на основе обычных мо- 
делей преобразования информации и приобретает вид, показанный в правой верхней 
части экрана на рис. 4 6). 

Топологическая картина ситуации, определяющей стремительный дрейф (рис. 5), 
принципиально отличается от обсуждаемой выше. Действительно, после восприятия 
удара происходит перестройка картины взаимодействия, расчет которой выполнен в 
рамках гипотезы квазистационарности. Судно, вовлеченное в стремительный дрейф, 
совершает эволюцию по пространственной траектории, проекция которой на плоскость 
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У, приведена на рис. 5. Здесь точками С1,С›,@з фиксируются положения ЦТ системы. 
Каждому такому состоянию системы соответствуют текущие значения характерных 
точек и областей топологической картины, типичной для надводного судна. ЦТ и ЦВ 
непрерывно изменяют свое положение. При этом формируется множество, отобража- 
ющее динамическую среду С7(60) и область изменения восстанавливающего момента 
М(6), которая трансформируется за счет многих неблагоприятных факторов вследствие 
сложной картины нелинейного взаимодействия, и в первую очередь — влияния воды на 
палубе от эффекта разрушения волны после удара. 


Рисунок 5 — Топологическая картина взаимодействия, отображающая движение 
судна в процессе эволюции под воздействием экстремальной волны 


Аналитическая интерпретация динамической картины взаимодействия в усло- 
виях стремительного дрейфа напоминает рассмотренную ранее картину для случая 
разрушающейся волны (рис. 4 6), но с другими характеристиками диаграммы остой- 
чивости М(0) и кренящего момента Мо. Однако принципиальные моменты взаимо- 
действия для рассматриваемой ситуации сохраняются: угол 0: здесь также характе- 
ризует неустойчивое положение равновесия при равенстве моментов М(0) и Мов 
точке 91, а дальнейшее увеличение крена 0 и опрокидывания происходят вследствие 
значительной величины кренящего момента и резкого снижения ординат диаграммы 
М(0). Проверка адекватности используемых математических моделей при анализе 
этой чрезвычайно сложной практической задачи ведется на основе специально 
разработанного подхода. Для этого был использован стенд, оборудованный в гидро- 
динамическом испытательном центре на основе функционирования виртуальной лабо- 
ратории интеллектуальной СКШ-системы [15]. 

Анализ альтернатив и окончательный выбор решения по управлению судном с 
соответствующей оценкой риска в каждом из разработанных сценариев взаимодей- 
ствия в рассматриваемой экстремальной ситуации ведутся на основе метода анализа 
иерархий [24] с использованием матрицы риска [25]. 

Работа выполняется в рамках проектов «Интеллектуальная система навигации 
и управления морским динамическим объектом в экстремальных условиях эксплу- 
атации», «Высокопроизводительный программный комплекс моделирования динамики 
корабля в экстремальных условиях эксплуатации» и «Интеллектуальные технологии 
поддержки процессов исследовательского проектирования судов и технических 
средств освоения океана» по направлениям «Создание интеллектуальных систем нави- 
гации и управления» и «Судостроение» в рамках программы «Научные и научно- 
педагогические кадры инновационной России» на 2009 — 2013 годы. 
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